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ABSTRACT 

Galaxy clusters provide us with important information about the cosmology of our universe. 
Observations of the X-ray radiation or of the SZ effect allow us to measure the density and 
temperature of the hot intergalactic medium between the galaxies in a cluster, which then allow 
us to calculate the total mass of the galaxy cluster. However, no simple connection between the 
density and the temperature profiles has been identified. Here we use controlled high-resolution 
numerical simulations to identify a relation between the density and temperature of the gas in 
equilibrated galaxy clusters. We demonstrate that the temperature-density relation is a real 
attractor, by showing that a wide range of equilibrated structures all move towards the attractor 
when perturbed and subsequently allowed to relax. For structures which have undergone sufficient 
perturbations for this connection to hold, one can therefore extract the mass profile directly from 
the X-ray intensity profile. 

Subject headings: galaxies: clusters: intracluster medium - methods: numerical 



1. Introduction 

Galaxy clusters are the largest known gravita 
tionally bound structures in the universe. Con- 
sisting almost entirely of dark matter and intr- 
acluster gas, they have proved to be very use- 
ful objects of study in cosmology. Cluster ob- 
servations have imposed const raints on the pa- 



rameters of the ACDM model Allen et all 12011 



IVikhlinin et al. 2009), and gravitational lensing 
reveals the galaxi es of the distant u niverse be- 
hind the clusters ( Kneib et al.ll2004r ). Observa- 



tions of the gas in galaxy clusters have enabled 
us to quantify the de n sity profile of dark matter 
dBuote fc Lewis! 12004 iPointecouteau et all [200l 
iHost fc Hansenll201ll ). However, both X-ray and 
SZ intensity observations suffer from the compli- 
cation that whereas it is easy to observe the inten- 
sity, which is a product of density and tempera- 
ture, then it is much more complicated to observe 
the t emperature and density separately (parazin 
Il986l ). Alternatively, if a simple relation between 
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the density and temperature would exist, one 
could extract the total mass profile directly from 
the easily observable X-ray or SZ intensities. 

Analytically one might try to search for such a 
relation by tu rning to the equation of hyd r ostati c 



19761) . 



equilibrium (jCavaliere fc Fusco-Femiano 
This equation relates the total mass of a clus- 
ter Mtot to the gas density p g and temperature T g 
of the intracluster gas. It is given by 



M tot (r) 
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where ks is the Boltzmann constant, fj, is the 
mean molecular mass, m u is the atomic mass 
unit, G is the gravitational constant and we define 
7 g = dln(p g )/dln(r) and r g = dln(T 9 )/dln(Y). 
Unfortunately, for a given M to t we can choose 
essentially any random p g (r) and still make the 
equation hold by finding a suitable T g (r). 

In this Letter we study the evolution of galaxy 
clusters using hydrodynamical simulations. We 
find that a wide range of initially equilibrated 
structures, when dynamically perturbed in or- 
der to mimic real perturbations happening during 
structure formation, all move towards a particular 
curve in the 2-dimensional phase space spanned 
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by the logarithmic derivatives of the density and 
temperature of the gas. 

2. Initial Conditions 

The dimensions of our structures are chosen to 
resemble those of a large galaxy cluster such as the 
Perseus Cluster. Each cluster has a scale radius of 
ro = lOOkpc and a total mass of 10 15 Mo of which 
10% is gas and 90% is dark matter. They are 
described by 1 million particles of which 10% rep- 
resent the gas and 90% represent the dark matter. 
The mass of the gas and the dark matter particles 
are thus the same, namely 10 9 M Q . 

To obtain a wide range of initial conditions, we 
create 20 different structures that vary in both 
gas and dark matter density profiles as well as 
dark matter velocity distributions. We use double- 
power laws to describe the density profiles of the 
dark matter and the gas. 
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Here a and 5 correspond to the inner and outer 
logarithmic slopes of the density profile respec- 
tively. The logarithmic slopes of the dark mat- 
ter density profiles are chosen in the range — 2 < 
7dm < in the inner region and — 6 < jdm < — 3.5 
in the outer region. For a given dark matter 
density profile we use the Eddington method to 
create isotropic dark matter velocity distributions 
with f3 = 0, or the Osipkov-Merrit method to cre- 
ate anisotropic dark matter velocity distributions 
with j3(r) = 1/(1 + (r a /r) 2 ), where r a = ro or 
r a = 2ro (Bin ney fc Tremaindliooi ). The loga- 
rithmic slopes of the gas density profiles are cho- 
sen in the range —1.5 < j g < in the inner re- 
gion and —5 < j g < —3.5 in the outer region. 
For a specific combination of gas and dark matter 
density profiles we use the equation of hydrostatic 
equilibrium to find the corresponding gas tempera- 
ture profiles. Note that 7 g and r g arc independent 
variables. It is therefore, in principle, possible to 
populate much of the space spanned by j g and r g . 

All of these initial conditions are presented in 
Figure [TJ To ensure stability these structures have 
all been evolved for 2 dynamical times, where the 
dynamical time td yn = 1.47Gyr has been com- 
puted at R = lOro. In all figures we only show 



Fig. 1. — The initial conditions in the "f g = 
dln( ( o s )/dln(r) and r g = dln(T g )/dln(r) space. 
The variations in gas and dark matter density pro- 
files as well as dark matter velocity distributions 
allow the structures to span a large portion of the 
T g -j g space. All the structures of this figure are 
stable on a timescale of 2 dynamical times when 
evolved in a hydrodynamical simulation. 



the region between 3 times the gravitational soft- 
ening and 10 times the scale radius. The values 
for T g and 7 g are the average of particles in loga- 
rithmically spaced radial bins. We emphasise that 
our initial conditions both include cool-core clus- 
ters with a central r g > 0, as well as non-cool core 
clusters with a central r g < 0. All these structures 
are stable on a timescale of 2 dynamical times. 



3. Numerical Code 

We use the parallel TreeSPH code G ADGET- 2 
to evolve the structures (|Springelll2005l ). GADGET- 
2 simulates fluid flows by means of smoothed par- 
ticle hydrodynamics and computes gravitational 
forces with a hierarchical tree algorithm. The time 
integration is based on a quasi-symplectic scheme, 
where long-range and short-range forces can be 
integrated with different time steps. Radiative 
cooling is implemented for a pr imordial mixture o f 
hydrogen and helium following iKatz et al.l (|1996l ) . 

The SPH properties are smoothed over the 
standard GADGET-2 kernel using 50 ±1 SPH par- 
ticles, while the gravitational forces are adjusted 
by the gravitational spline kernel using a soften- 
ing length of 6.0kpc. To protect against too large 
time steps for particles with very small acceler- 
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ations, we require that the maximum time step 
is two percent of the dynamical time of the sys- 
tem. We use an artificial viscosity parameter of 
oivisc = 0.8 and a courant factor of a cour = 0.1. 

For modelling star formation and the asso- 
ciated heating by supernovae (SN) we follow 
the sub-r esolution multiphase I SM m odel devel- 
oped by ISpringel k Hernquistl (|2003l ). In this 
model, a thermal instability is assumed to oper- 
ate above a critical density threshold pth, pro- 
ducing a two-phase medium which consists of 
cold clouds embedded in a tenuous gas at pres- 
sure equilibrium. Stars are formed from the cold 
clouds on a times cale chosen to match observations 
(jKennicuttl Il998l ) and short-lived stars supply an 
energy of 10 51 ergs to the surrounding gas as SN. 
We adopt the standard parameters for the multi- 
phase model i n order to match the K e nnicu tt Law 
as specified in ISpringel k Hernquistl (|2003l ). The 
star formation timescale is set to = 2.1Gyr, 
the cloud evaporation parameter to Ao = 1000 
and the SN "temperature" to T SN = 10 8 K. More 
information on the GAD GET version adopt ed in 
this work can be found in lMoster et al.l (|201lh . All 
simulations are carried out in a non-cosmological 
Newtonian box. 

4. Perturbations 

The cosmological structures we observe today 
have a long history of hierarchical structure forma- 
tion involving gravitational collapse and mergers. 
Both of these are violent relaxation processes that 
change the gravitational potential and therefore 
also the energy of the particles according to 



dE 
dt 
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for spherical systems, where E is the energy, $ 
i s the gravitational potential and t is the time 
l|Lvnden-Bei]||l967t iBinnev k Tremainell2008l) 

To mimic these processes we will expose all our 
initially equilibrated structures to controlled per- 
turbations. We will require that these perturba- 
tions are both continuous and spherical. One ex- 
ample of such perturbations is to let the value of 
the g ravitational constant G vary in the N-body 
code (jSparre k Hansenll2012h . The perturbations 
are spherical since the structures themselves are 
spherical. They are also continuous since they af- 



fect the accelerations and not the velocities of the 
particles. 

Initially, we change the value of G by 10%. Set- 
ting the gravitational constant equal to 1.1G we 
run the structures for one dynamical time. This 
will cause the structures to contract, since a higher 
value of G results in a deeper potential. We then 
set the gravitational constant equal to 0.9G and 
run the structures for another dynamical time. 
This will cause the structures to expand since a 
lower value of G results in a more shallow poten- 
tial. All in all we submit each structure to a series 
of 20 perturbations, alternating between a larger 
and a smaller value of the gravitational constant as 
described. We check that no unwanted oscillatory 
behaviour results. We then proceed with ten per- 
turbations where the value of G is changed by 5%, 
and then ten perturbations where the value of G is 
changed by only 1%. Finally, we let the structures 
run for an additional 10 dynamical times with a 
normal value of G to make sure they are in equi- 
librium. 



-2 r 




Fig. 2. — The same phase space as in Figure [T] 
is now populated by the final equilibrated profiles 
of the perturbed structures. Note that the initial 
conditions (light grey symbols) have approached 
the 1-dimensional line, the attractor, from both 
sides. Also shown is a line of the form p g jT^ 2 ~ 
r~ 2 (dashed), and a linear guide-the-eye line of the 
form 7 9 = 0.83r g — 1.72 in the region —0.7 < r g < 



0.6 (solid). 
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5. The Attractor Solution 

The results of the perturbations are presented 
in Figure[5J from which it is clear that all the struc- 
tures end up along the same curve. It is important 
to emphasise that we had initial conditions both 
above and below this attractor solution. The par- 
ticles in the outer region (small r) still tend to 
follow the initial conditions. This is partly be- 
cause they have not been perturbed as much due 
to their larger dynamical time. 

For the final structures the formed stars dom- 
inate the total mass within r = 20kpc, and they 
contribute significantly to the collisionless compo- 
nent within r — lOOkpc. To determine how much 
of an effect the star formation has had on the at- 
tractor, we first perform a simulation in which 
the density threshold for star formation, p t h, is 
ten times larger than normal. This has the ef- 
fect of suppressing the star formation. However, 
the structure still ends up along the attractor sug- 
gesting that the star formation itself does not have 
much of an effect. 

We also perform a run of only 20 perturbations 
of 10%. After that we let the structure relax for 
30 dynamical times. Again, the structure does not 
stray from the attractor at all, showing that it is 
not affected by the formation of stars even over 
such a long period of time. 

To quantify how much the attractor solution 
stands out compared to the initial conditions 
we introduce the structure-to-structure dispersion 
c| 2 s( r s); which describes the spread among a 
group of structures as a function of r g . This value 
is expected to be significantly larger for the initial 
conditions than for the perturbed structures, be- 
cause the different structures make clearly distinct 
lines in Figure [TJ One can define 



a tot( T g) = a S2s( T g) + °fni 



(r g ), 



(4) 



where o\ ot {rg) is the total spread of data points, 
and <jf nt {T g ) is the average internal spread of the 
individual structures. Perfectly smooth structures 
with infinitely many particles would thus have 
a int( T g) 1=3 0- Binning the data points with re- 
spect to T g we consider a single bin i with N data 
points. Since the data points are not distributed 
symmetrically around an average value of j g , we 
define the median value of 7 g as the one for which 
50% of the data points lie to either side. Simi- 



larly, counting 15.8% from the left (one standard 
deviation) we find 

cr_(rg) = (50— 15.8 percentile) -7* (50 percentile)!, 

(5) 

and counting 15.8% from the right we find 
cr+(Tg) = |7*(50+15.8 percentile) -7* (50 percentile)). 

We now define 



a tot\ T q) = o 



(0) 
(7) 



and its error 

S°L(tD = \{{oArl)-^ot{rl)) 2 + {a + {rl)-aUrl)) 2 ) 
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For a single structure with M radial bins we define 



A lg =ll-% 
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as the distance between a point 7^ and its inter- 
polated value 



V+ 1 



v'- 1 
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Here j refers to a given radial bin. Again, binning 
in T g we define for a single bin i with K data points 
the value of c? ni (r*) as the median value of (A7 ff ) 2 . 
For the distribution we again count 15.8% from 
left and right and take the average to to find the 
error Sa^ nt (Tg). We do not include any points for 
which 7j < —3.5. 

The spread between the initial and final struc- 
tures are presented in Figure |3l Since the initial 
conditions span a large portion of the 7 g — r g space 
they consequently have a relatively large cr|. 2S , 
while the final structures, which generally lie along 
the same curve, have a distinctly smaller c^s- 



6. Discussion 
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In Figure[5]we plot a line of the form p g /T g 
r~~ a , corresponding to the entropy following a 
power-law in radius, where a is an unknown con- 
stant. This implies 7 9 = —a + 3/2 r g . This possi- 
bility has bee n discussed intensively in the litera- 
ture (see e.g. iMorandi k, Ettori|[2007t iPratt et al 
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Fig. 3. — The structure-to-structure dispersion, 
describing the off-set between a group of different 
structures. The initial conditions (blue) cover a 
large fraction of the phase space (T g ,^ g ), whereas 
the final (green) structures end essentially on top 
of each other, resulting in a very small <J%2S- 



20101 Dubois et al. 2011 . for references). In the 
figure we use a = 2. Our resulting attractor 
clearly differs from this simple behaviour, both in 
the inner and the outer region. 

We also include a linear guide-the-eye line of 
the form j g = 0.83t 9 — 1.72. Although this is 
obvious a simplified approximation it provides an 
acceptable fit in the range —0.7 < r g < 0.6. 

We note that our limited resolution does not al- 
low us to resolve the inner region where 7 S > — 1 . 
That region is particularly interesting when one 
is considering cool-core clusters where r g may be 
large and positive. Initially the rapid cooling in 
our simulations forms a large central stellar con- 
densation, which acts like a (unrealistically large) 
central blue BCG, a problem which possibly could 
be reduced by including AGN feedback. Subse- 
quently, our perturbations are likely inducing adi- 
abatic compression shocks, which helps (together 
with the SN feedback) preventing the structures 
from undergoing further rapid cooling, in fact, 
we do not observe any late time cooling catastro- 
phe despite our simulation effectively being longer 
than a Hubble time. This may also be connected 
to the fact that the quickly formed central BCG 
may act like a stabilizing structure which helps 
preventing ra pid over cooling f rom destabilizing 
the structure (|Li fe Brvanll2011l ). 



This simple connection allows one to write 
a simplified hydrostatic equilibrium equation, 
M to t = fi(T g ), where fi(T g ) is a known function 
only of the gas temperature, with no reference 
to the gas density. Alternatively, if we knew the 
absolute value of the temperature at one point, 
then one could write this as M tot = f2(p g ), where 
fiiPij) is a known function only of the gas density, 
or even better, M to t = /a CO, where /s(/) is a 
known function only of the gas intensity, which is 
easily observable even at high redshift. 

It has long been observed that the characteris- 
tic radius for the temperature (at T max ) roughly 
corresponds to the scale radius fo r the gas den- 
sity (see e.g. IVikhlinin et al.l 120061 ) . We now see 
how this relates to the attractor presented here: 
The radius corresponding to r g = should have 
—2 < "f g < —1.5, which is close to the scale radius 
of the density profile 

A connection has previously been suggested be- 
tween the gas tempe rature profile and t h e dar k 
matter dispersion (jHansen fc Piffarettil 120071 ). 
which is not only a fair approximation for galax 



wnicn is not only a lair approximation lor galaxy 
clusters jFTederiksen et al.ll2009HHost et al.lbooa 
ZuHonell201ll ). but also holds well f or the hot gas 



in ga laxies as well as near AGN ([Hansen et al 



201 if ). We find that the corresponding structure- 
to-structure dispersion for this gas-DM tempera- 
ture connection is larger, and we therefore con- 
clude that the gas-density connection presented in 
this paper is likely more fundamental. 

When comparing our finding here with the re- 
cent discovery of an attractor for dark matter 
structures, a remarkable suggestion presents it- 
self. The dark matter attractor is a 1-dimensional 
line in a 3-dimensional space spanned by 7dm, Tdm 
and Pdm, where the first two are the logarith- 
mic derivatives of the dark matter density and 
radial velocity dispersion (of, which acts like a 
temperature p er mass), and 0dm. is the veloc 



ity anisotropy ( Hansen. Juncher fc Sparre 2010l) . 
The Jeans equation, which is the dark matter 
equivalent of the hydrostatic equilibrium, can be 
written as 

2 

M tot (r) = -r^ (7dm(r) + r dm {r) + 2f3 dm (r)) . 

(11) 

We therefore see that the dark matter attractor is 
between elements in the parenthesis in the Jeans 
equation, eq. (|TT]) , and the gas attractor presented 
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in this paper is between the elements in the paren- 
thesis in the equation for hydrostatic equilibrium, 
eq. (HJ). We therefore speculate, that if one will 
manage to derive the gas attractor we present 
here, then one may possibly use that understand- 
ing to derive the attractor for the dark matter as 
well. 

It is important to emphasize that the attractor 
for the gas is not just an reaction to the changes 
in the dark matter as it adjusts to its own attrac- 
tor. Two of our initial structures start right on 
the dark matter attractor but not on the gas at- 
tractor. When they are perturbed they move to 
the gas attractor while staying on the dark matter 
attractor. 

In this paper we have only investigated the ef- 
fect of the particular G-perturbation. We note 
that the resulting dark matter density profiles are 
different, and the gas attractor is thus not a trivial 
result of the perturbation driving the dark mat- 
ter structures towards a universal density profile. 
It has been shown for the dark matter attractor 
that a wide ran ge of perturbations all lead to the 
same attractor (Hansen Juncher k Sparrd 2010t 
Sparre k Hansen! l2012t iBarber et al.ll2012i sug- 
gesting that our found gas attractor is likely not a 
result of the particular perturbation chosen here. 
In the future we intend to study which effects more 
realistic perturbations, such as minor or major 
mergers, will have on the gas attractor. We also 
note that the simulations presented here both have 
gas and dark matter, however, when analysing the 
3-dimensional dark matter phase space, we find 
that the dark matter attractor is slightly shifted 
compared to the cases including only dark matter. 

7. Conclusions 

Using high-resolution numerical simulations of 
large galaxy clusters, including radiative cooling, 
star-formation and feedback, we identify an at- 
tractor for the hot gas component. This attractor 
gives a very simple connection between the tem- 
perature and the density of the gas. 

We find the attractor to hold for all structures 
which have been exposed to significant perturba- 
tions and subsequently allowed to relax. Since cos- 
mological structures experience perturbations dur- 
ing mergers, we speculate that the central part of 
galaxy clusters, which are fully equilibrated, may 



follow the attractor. It will be very interesting 
in the near future to test this suggestion on real 
X-ray observations. 
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Research Foundation. 

REFERENCES 

Allen, S. W., Evrard, A. E., k Mantz, A. B. 2011, 
ARA&A, 49, 409 

Barber, J. et al. 2012, to appear 

Binney, J., k Tremaine, S. 2008, Galactic Dynam- 
ics: Second Edition. Princeton University Press 

Buote, D. A., k Lewis, A. D. 2004, ApJ, 604, 116 

Cavaliere, A., k Fusco-Femiano, R. 1976, A&A, 
49, 137 

Dubois, Y., Devriendt, J., Teyssier, R., & Slyz, A. 
2011, MNRAS, 417, 1853 

Frederiksen, T. F., Hansen, S. H., Host, O., & 
Roncadelli, M. 2009, ApJ, 700, 1603 

Hansen, S. H., & Piffaretti, R. 2007, A&A, 476, 
L37 

Hansen, S. H., Juncher, D., & Sparre, M. 2010, 
ApJ, 718, L68 

Hansen, S. H., Maccio, A. V., Romano-Diaz, E., 
et al. 2011, ApJ, 734, 62 

Host, O., Hansen, S. H., Piffaretti, R., et al. 2009, 
ApJ, 690, 358 

Host, O., & Hansen, S. H. 2011, ApJ, 736, 52 

Katz, N., Weinberg, D. H., & Hernquist, L. 1996, 
ApJS, 105, 19 

Kcnnicutt, R. C, Jr. 1998, ApJ, 498, 541 

Kneib, J.-P, Ellis, R. S., Santos, M. R., & 
Richard, J. 2004, ApJ, 607, 697 

Li, Y., k Bryan, G. L. 2011 



6 



Lynden-Bell, D. 1967, MNRAS, 136, 101 

Morandi, A., & Ettori, S. 2007, MNRAS, 380, 
1521 

Moster, B. P., Maccio, A. V., Somerville, R. S., 
Naab, T., & Cox, T. J. 2011, MNRAS, 415, 
3750 

Pointecouteau, E., Arnaud, M., & Pratt, G. W. 

2005, A&A, 435, 1 

Pratt, G. W., Arnaud, M., Piffaretti, R., ct al. 
2010, A&A, 511, A85 

Sarazin, C. L. 1986, Reviews of Modern Physics, 
58, 1 

Sparre, M., & Hansen, S. H. 2012, to appear 

Springel, V., & Hernquist, L. 2003, MNRAS, 339, 
289 

Springel, V. 2005, MNRAS, 364, 1105 

Vikhlinin, A., Kravtsov, A., Forman, W., et al. 

2006, ApJ, 640, 691 

Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et 
al. 2009, ApJ, 692, 1060 

ZuHone, J. A. 2011, ApJ, 728, 54 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



7 



